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Using a self-consistent Bogoliubov-de Gennes approach, we theoretically study the proximity-induced den¬ 
sity of states (DOS) in clean SFF spin-valves with noncollinear exchange fields. Our results clearly demon¬ 
strate a direct correlation between the presence of a zero energy peak (ZEP) in the DOS spectrum and the 
persistence of spin-1 triplet pair correlations. By systematically varying the geometrical and material param¬ 
eters governing the spin-valve, we point out to experimentally optimal system configurations where the ZEPs 
are most pronounced, and which can be effectively probed via scanning tunneling microscopy. We complement 
these findings in the ballistic regime by employing the Usadel formalism in the full proximity limit to investigate 
their diffusive SFF counterparts. We determine the optimal normalized ferromagnetic layer thicknesses which 
result in the largest ZEPs. Our results can serve as guidelines in designing samples for future experiments. 

PACS numbers: 74.50.+r, 74.25.Ha, 74.78.Na, 74.50.+r, 74.45,+c, 74.78.FK, 72.80.Vp, 68.65.Pq, 81.05.ue 


I. INTRODUCTION 

The interplay of ferromagnetism and superconductivity in 
hybrid superconductor ( S ) ferromagnet (F) structures (S/F 
structures) constitutes a controllable system in which to study 
fundamental physics, including prominently that of compet¬ 
ing multiple broken symmetries. 1 ' 2 The proximity of a con¬ 
ventional s-wave superconductor with non-aligned ferromag¬ 
netic layers, or a textured ferromagnet, induces both spin- 
singlet and odd frequency 3 (or equivalently odd-time 4 ) spin- 
triplet correlations with 0 and ±1 spin projections along a 
spin quantization axis. These triplet pairs stem from bro¬ 
ken time reversal and translations 1,2 symmetries. This kind 
of spin-triplet pairings originally suggested as a possible pair¬ 
ing mechanism in 4 He, 3 has reportedly been observed in in- 
termetallic compounds such as Sr 2 Ru 04. 5,6 SF heterostruc¬ 
tures are particularly simple and feasible experimental sys¬ 
tems which allow for direct studies of the intrinsic behavior of 
differing superconducting pairings. Unlike the opposite-spin 
correlations, spin-1 pairing correlations are rather insensitive 
to the pair-breaking effects of ferromagnetic exchange split¬ 
ting, and hence to the thickness of the magnetic layers, tem¬ 
perature, and magnetic scattering impurities. The amplitudes 
of the opposite-spin correlations pervading the adjacent fer¬ 
romagnet, undergo damped oscillations as a function of po¬ 
sition which reveals itself in 0-7T transitions of the supercur¬ 
rent. 17 10 Since about a decade ago, several proposals have 
been put forth to achieve attainable and practical platforms 
that isolate and utilize the proximity-induced 14 superconduct¬ 
ing triplet correlations in SF hybrids. 12 

The signatures of the proximity-induced electronic den¬ 
sity of states (DOS) in the F layers of these hybrid struc¬ 
tures can reveal the existence and type of superconducting 
correlations in the region. 1,-1 3 One promising prospect for un¬ 
ambiguously detecting triplet correlations experimentally in¬ 
volves tunneling spectroscopy experiments which can probe 
the local single particle spectra encompassing the proximity- 
induced DOS. 1112,14 29 Nonetheless, competing effects can 
make analysis of the results of such a ‘direct’ probe of spin- 


triplet superconducting correlations problematic. The DOS in 
SNS junctions and SFS heterostmctures where the magne¬ 
tization pattern of the F layer can be either uniform or tex¬ 
tured (including domain wall and nonuniform textures, such 
as the spiral magnetic structure of Holmium) has been ex¬ 
tensively studied. 11,12,20-22,30,31 It was found that the DOS in 
a normal metal sandwiched between two s-wave supercon¬ 
ducting banks shows a minigap which closes by simply tun¬ 
ing the superconducting phase differences up to the value of 
n 11,12,21,22 

In contrast, the DOS can exhibit anomalous behav¬ 
ior in inhomogeneous magnetic layers. Namely, upon modu¬ 
lating the superconducting phase difference 11,12 a peak arises 
at zero energy, at the center of what was a minigap. It was 
also shown that the zero energy peak (ZEP) in the DOS for 
a simple textured SFS junction can be maximized at a 7r 
bias. 11,12 The minigap-to-peak behavior of the DOS at zero 
energy is an important signature of the emergence of triplet 
correlations. 11,12 Recently it was theoretically proposed that 
the minigap-to-peak phenomenon be leveraged for functional¬ 
ity in device platforms such as SQUIDs, to enhance their per¬ 
formance and as ultrasensitive switching devices, including a 
singlet-triplet superconducting quantum magnetometer. 23 

An important spectroscopic tool for investigating proxim¬ 
ity effects on an atomic scale with sub-meV energy resolu¬ 
tion is the scanning tunneling microscope (STM). As shown 
in Fig. 1, an SFF spin valve structure can be probed exper¬ 
imentally by positioning a nonmagnetic STM tip at the edge 
of the sample to measure the tunneling current (I) and volt¬ 
age (V) characteristics. This technique yields a direct probe 
of the available electronic states with energy eV near the tip. 
Therefore, the differential conductance dI(V)/dV over the 
energy range of interest is proportional to the local DOS. Nu¬ 
merous experiments have reported signatures of the energy 
spectra in this manner. 14,15,17-19,24-26,28 When ferromagnetic 
elements are present, the superconducting proximity-induced 
DOS reveals a number of peculiarities due to the additional 
spin degree of freedom that arises from the magnetic layers. 
However, the experimental signatures of the odd-frequency 
spin-triplet correlations can be washed out by more dominate 
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singlet correlations. When the exchange splitting h of the 
magnetic layers is large (~£f, i.e. close to the half metal¬ 
lic limit), the characteristic length scale £p that describes the 
propagation length of opposite-spin pairs in the ferromagnets 
is extremely small. These types of proximity-induced correla¬ 
tions can thus only be experimentally observed in weak mag¬ 
netic alloys h <C £p (such as Cu r Ni y ) or thin F layers so that 
is sufficiently large to allow the opposite-spin super¬ 
conducting correlations to propagate in the ferromagnet with¬ 
out being completely suppressed. 14,16 - 27 Since spin-1 triplet 
pairs are not destroyed by the ferromagnetic exchange field 
in strong magnets, there should exist certain system parame¬ 
ters, e.g., ferromagnet widths and exchange fields, that result 
in regions whereby equal-spin pairs are the only pair correla¬ 
tions present. This scenario was explored in a S'/Ho bilayer 29 , 
where phase-periodic conductance oscillations were observed 
in Ho wires connected to an ordinary s-wave superconductor. 
This behavior was qualitatively explained in terms of the long- 
range penetration of proximity-induced spin-1 triplet pairings 
due to the helical structure of the magnetization. 2 In prac¬ 
tice however, simpler structures involving SF hybrids with 
uniform exchange fields are often preferable from both an ex¬ 
perimental and theoretical perspective. I4 ’ 16,17,24-2y Therefore, 
the primary aim of this work is the determination of experi¬ 
mentally optimal parameters for probing odd-frequency spin- 
1 triplet correlations with DOS signatures in nanoscale SFF 
spin valves. 

Nearly all of the past theoretical works on SFF struc¬ 
tures have considered the diffusive case, 33-35 where impurities 
strongly scatter the quasiparticles. The clean regime has been 
studied, using a self-consistent solution of the Bogoliubov- 
de Gennes (BdG) 39 equations, in Ref. 37. That work, how¬ 
ever, focused largely on the transition temperature oscilla¬ 
tions. The results for these oscillations were found 38 to agree 
with experiment and to be consistent with other experimen¬ 
tally established 39 42 results. In the present work, we use 
the same general methods used there to study a simple SFF 
structure with noncollinear exchange fields in the ballistic 
regime, but we focus on a very different quantity which is 
readily accessible experimentally, namely the local DOS and 
its detailed low-energy structure. We strongly emphasize the 
relation between the ZEP and the triplet pairing amplitudes. In 
particular, given the assertion ' 8 that variations in the transition 
temperature in these valve stmctures are quantitatively related 
to the average triplet pair amplitudes in the outer F layer, we 
will search for, (and, as will be seen, find) correlations be¬ 
tween the ZEP and these averages. This BdG study is com¬ 
plemented with a briefer investigation of the corresponding 
diffusive case. By considering both regimes, we will be able 
to provide some general guidelines for future experiments. 

The structure we study is schematically depicted in Fig. 1, 
where the STM tip is positioned at the outermost F layer, 
near the vacuum boundary. In the ballistic regime, we employ 
the full microscopic BdG equations within a self-consistent 
framework. From the solutions, we calculate the local DOS 
over a broad range of experimentally relevant parameters, and 
study its behavior at low energies. For the diffusive regime, 
we make use of the quasiclassical Usadel 43 approach to study 


the diffusive SFF counterparts in the full proximity limit. 
Our systematic investigations thus provide a comprehensive 
guide into such spin valves. Utilizing experimentally realis¬ 
tic parameters, we determine favorable thicknesses for the F 
layers to induce maximal ZEPs, which occurs when the pop¬ 
ulation of triplet correlations in the outer layer, dominates the 
singlets. 

The paper is organized as follows. In Sec. II we outline 
the theoretical approaches used. In Sec. Ill, we present our 
results in two subsections, pertaining to the ballistic and dif¬ 
fusive regimes. In the ballistic case, we study the local DOS 
for differing exchange field misalignments, exchange field in¬ 
tensities, and interface scattering strengths. We also investi¬ 
gate the singlet and triplet pairing correlations for similar pa¬ 
rameters to determine how ZEPs in the DOS correlate with 
the triplet correlations. In the diffusive case, we present two- 
dimensional maps of the ZEP at different exchange field mis¬ 
alignments, and SF interface opacity. Finally, we summarize 
with concluding remarks in Sec. IV. 


II. METHODS AND THEORETICAL TECHNIQUES 


In this section, we first discuss the theoretical framework 
used to study clean samples. We then outline the Usadel tech¬ 
nique in the full proximity regime, which properly describes 
dirty samples. 



FIG. 1. (Color online) Schematic of the SFF spin-valve structure. 
The ferromagnetic layers have uniform exchange fields located in 
the yz plane. The exchange field of each layer is defined by h 1,2 = 
ho(0, sin /3 i,2, cos A, 2 ) in which /3 1,2 are the angle of the exchange 
fields with respect to the 2 direction. The ferromagnets (Fi, FS) 
and superconductor (S) are stacked in the x direction with thickness 
dp 1 , dF 2 , and ds, respectively. The STM tip is located at edge of 
the SFF spin-valve. 
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A. Microscopic approach: Bogoliubov-de Gennes equation 


For the ballistic regime, we use the microscopic BdG equa¬ 
tions to study SFF spin valve nanostructures. We solve these 
equations in a fully self-consistent 37,44 manner. A schematic 
of the spin valve configuration is depicted in Fig. 1. The gen¬ 
eral spin-dependent BdG equations for the quasiparticle ener¬ 
gies, e n , and quasiparticle amplitudes, u na , v na is written: 
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where the pair potential, A(x), is calculated self-consistently 
as explained below. This quasi one-dimensional system is de¬ 
scribed by the single particle Hamiltonian 'Ho(x) as, 

n 0 (x) = ^ {-dl + k 2 y + k 2 z ) -e f + u{x), (2) 

where Ep is the Fermi energy, and U(x) is the spin- 
independent interface scattering potential which we take to 
be of the form U(x) = F[[5{x — dpi) + 5(x — dpi — dp 2 )]. 
The in-plane wavevector components, k y and k z , arise from 
the translational invariance in the y and z directions. The sys¬ 
tem is finite in the x direction, with widths of each F and S 
layer shown in the schematic. Our method permits arbitrary 
orientations and magnitudes of the magnetic exchange fields, 
hi (i = 1,2), in each of the ferromagnet regions. Specifically, 
we fix the exchange field in F 2 to be aligned in the z direction, 
while in F-\ , its orientation is described by the angle /?i: 

^ _ fdi = d 0 (0,sin/3i,cos/3i), in 
\h 2 = h 0 z, in F 2 , 

where we consider the experimentally appropriate situation of 
an in-plane Stoner-type exchange field interaction. 

The spin-splitting effects of the exchange field coupled with 
the pairing interaction in the S regions, results in a nontrivial 
spatial dependence of the pair potential A(x). In general, it 
is necessary to calculate the pair potential in a self consistent 
manner by an appropriate sum over states: 

A (z) = + Uni(x)v„ t {x)\ tanh , 

n 

(4) 


where g(x) is the attractive interaction that exists solely inside 
the superconducting region and the sum is restricted to those 
quantum states with positive energies below an energy cutoff, 

Wp>. 

We now discuss the appropriate quantities that characterize 
the induced triplet correlations. We define 4,43 the following 


triplet pair amplitude functions in terms of the field operators 
in the Heisenberg picture, 

fo(x,t) = ^ [(^(MXOgO)) + (ipi(x,t)ip t {x, 0))], 

(5a) 

fi{x,t) = \ (x,0)) - (VtOMXOg 0 ))], 

(5b) 

where t is the relative time. With the quantization axis aligned 
along the z direction, the time-dependent triplet amplitudes, 
fo(x,t) and f\(x,t), can be written in terms of the quasipar¬ 
ticle amplitudes: 4,43 

fo(x,t) = ( fn l ( x ) - fn( x )) C n(t), (6) 

n 

fl(x,t) = \Y^ + f^( X )) C n(t), (7) 

n 

where we define f° a {x) = u na (x)Vn a i {x), and the time fac¬ 
tor ( n {t) is written, 

( n (t) = cos(e n t) - i sin(e„f) tanh ■ (8) 

Experimentally accessible information regarding the quasi¬ 
particle spectra is contained in the local density of one particle 
excitations in the system. This includes the zero-energy sig¬ 
natures in the density of states (DOS), which present a pos¬ 
sible experimental avenue in which to detect the emergence 
of equal-spin triplet correlations within the outer ferromag¬ 
net. The total DOS, N{x, e), is the sum N^(x. e) + N±(x, e), 
involving the spin-resolved local density of states (DOS), N a , 
which are written, 

N a{ x,e) = -^{[<(a:)] 2 / , (e-e„) + K(a:)] 2 / , (e+e rl )}, 

n 

(9) 

where er denotes the spin (=1\4-X an d f{c) = df/de is the 
derivative of the Fermi function. 


B. Quasiclassical approach: Usadel equation 


When the system contains a strong impurity concentration, 
then for sufficiently small energy scales, the superconducting 
correlations are governed by the Usadel equation. Following 
Ref. 12, the Usadel equation 1 compactly reads: 


V d,G{r,e)[d,G(r,e)] 
cp 3 + diag [h(r ) • er, (h(r) 


>r-i 


) G(r , e) 


= 0 ,( 10 ) 


in which p 3 and er = [a x , a v , a z ) are 4 x 4 and 2x2 Pauli ma¬ 
trices, respectively, and V represents the diffusive constant of 
the magnetic region. The quasiclassical approach employed 
in this section supports ferromagnets with arbitrary exchange 
field directions; h(r) = (h x (r), h v (r) , h z (r)). In Eq. (10), 








4 


G represents the total Green’s function which is made of Ad¬ 
vanced (A), Retarded (R), and Keldysh (K) blocks. There¬ 
fore, the total Green’s function can be expressed by: 

G{r,e) = (^' Q G gA ) , G R (r,e) = ^ . 

( 11 ) 


In the presence of ferromagnetism, the components of ad¬ 
vanced block, G A (r ), of total Green’s function G can be writ¬ 
ten as: 


^-(&!*)■ S(r ' e) =fe»u)- (12) 

In this paper, however, we assume stationary conditions for 
our systems under consideration, and hence the three blocks 
comprising the total Green’s function are related to each other 
in the following way: G A (r,e) = — [p3G R (r, e)p%\ , and 
G A (r,e) = tanh(/3e) [G R (r, e) — G j4 (r,e)], where j3 = 
k B T/2. 

The SF interface controls the proximity effect. There¬ 
fore, appropriate boundary conditions should be considered 
to properly model the system. In our work, we consider 
the Kupriyanov-Lukichev boundary conditions at the SF 
interface 46 which controls the induced proximity correlations 
using a parameter £ as the barrier resistance: 

(G(r , e)dG(r , e) = [G BC s(0, e), G(r, e)]. (13) 

The solution for a bulk even-frequency s-wave superconduc¬ 
tor G R CS reads, 32 


/ lcoshi?(e) icr y sinfn?(e) \ 

^ ia y sinht?(e) — lcosh$(e) J’ 


(14) 


where i}(e) = arctanh(| A | /e). 

The system local density of states, A f(r, e), can be ex¬ 
pressed by the following equation: 




•A/o 

T 


Re 


Tr{G(r, e)} 


(15) 


in which Ao is the density of state normal state. 


III. RESULTS AND DISCUSSION 

In this section, we describe our results. We start with those 
for a ballistic SFF structure and then present the predictions 
of Usadel formalism for diffusive samples. 


A. Ballistic Regime 

In this subsection we present the self-consistent results for 
the ballistic regime. The numerical method used here to it¬ 
eratively solve in a self consistent way Eqs. (1) and (4) has 
been extensively described elsewhere, 37 - 44 and details need 
not be repeated here. In the calculations, the temperature 


T is held constant at T = 0.05T C , where T c is the transi¬ 
tion temperature of a pure bulk S sample. All length scales 
are normalized by the Fermi wavevector, so that the coordi¬ 
nate x is written X = kpx, and the F\ and F 2 widths are 
written Dpt = kpdpi, for i = 1,2. The ferromagnet F 2 
and superconductor are set to fixed values, corresponding to 
Dp 2 = 400, and Dg = 600, respectively. We also assume a 
coherence length corresponding to kpt; o = 100. One of our 
main objectives in this paper is to study the triplet correla¬ 
tions, which are odd in time. 4 To accomplish this, we employ 
the expressions in Eqs. (6) and (7), which describe the spa¬ 
tial and temporal behavior of the triplet amplitudes. At t = 0 
the triplet correlations vanish because of the Pauli exclusion 
principle. At finite t, the triplet correlations generated near 
the S/F interface tend to increase in amplitude and spread 
throughout the structure. We normalize the time t according 
to r = LOpt, and we set it to a representative 37 value of r = 4. 
We can then study the behavior of the triplet amplitudes /o and 
/i throughout the junction. To explore the proximity induced 
signatures in the single-particle states, which is the main pur¬ 
pose of this work, we then present a systematic investigation 
of the experimentally relevant local DOS. All DOS results 
presented are local values taken at a fixed position near the 
edge of the sample in the 7*2 region. We characterize inter¬ 
face scattering, when present, by delta functions of strength 
H, which we write in terms of the dimensionless parameter 
H b = H/vf. Finally, we use natural units, e.g., h = k B = 1 
throughout. 


Triplet and singlet pair correlations 

Here we present results for both the triplet and singlet cor¬ 
relations, calculated using Eqs. (6)-(7). For the cases shown 
below, the absolute value of the singlet and triplet complex 
quantities are averaged over the region of interest, which in 
this case is the experimentally probed F 2 region. An impor¬ 
tant reason for focusing on those spatially averaged (over the 
outer magnet) quantities, rather than the spatial profiles dis¬ 
cussed in Ref. 37, is that it was experimentally shown’ 3 that 
these triplet averages perfectly anticorrelate with the transi¬ 
tion temperatures, i.e. the spin valve effect. We also normal¬ 
ize all pair correlations to the value of the bulk singlet pair 
amplitude. We begin by showing, in Fig. 2, the spatially aver¬ 
aged absolute value of the complex triplet amplitudes |/o,avg| 
(with spin projection m = 0), and |/i, a vg| (with spin pro¬ 
jection m = ±1) along with the singlet |/ 3 )avg |, (note that 
fz(x) = A(x)/g(x)) as functions of Dpi. Each row of pan¬ 
els corresponds to a different exchange field value: from top to 
bottom rows, we have h/sF = 0.5,0.1, and 0.05. Examining 
the opposite spin correlations, fo and damped oscillatory 
behavior with Dpi is evident: this is related to the spatial os¬ 
cillation of the Cooper pairs due to their acquiring a center of 
mass momentum when entering the magnet. 47 Therefore, the 
wavelength of these oscillations varies inversely with the ex¬ 
change field in F 2 (this is why the Dpi range for the weaker 
exchange fields is extended). Quantum interference effects 
generate peaks in /o and that occur approximately when 
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FIG. 2. (Color online) The absolute value of the normalized triplet 
and singlet pair correlations, averaged over the F 2 region, as a func¬ 
tion of Dfi. The exchange field strengths are (from top to bottom): 
h/eF = 0.5,0.1,0.05. The relative exchange field orientations are 
orthogonal, with 0 i = tt/2, and 02 = 0. 


^fi/£f = nn, (i.e. Dpi = nn{h/£p)~ l ). In the ballistic 
regime, the length scale that characterizes the damped oscil¬ 
lations is = vf/ (2 h), where 'up is the Fermi velocity. The 
equal-spin amplitudes / 1 , are seen to behave oppositely, with 
a phase offset of approximately n/2. Their magnitude de¬ 
clines more rapidly with Dpi, compared to the behaviors of 
/o and / 3 . This is consistent with /1 triplet generation being 
optimal for highly asymmetric ferromagnetic layer widths. 32 
It is notable that the periodic occurrence of peaks in f\ when 
varying Dpi, evolves into a single maximum as h is reduced 
further. 

One of the strengths of the microscopic BdG formalism 
is having the ability to properly include the full microscopic 
range of length and energy scales inherent to the problem. 
This includes the exchange energy h, which in our BdG 
framework can span the limits from a nonmagnetic normal 


FIG. 3. (Color online) The spatially averaged (in F 2 ) normalized 
triplet and singlet pair correlations as a function of Ii/sf- As in 
Fig. 2, the magnitude of each quantity is taken and averaged over 
the F 2 region. Each row of panels corresponds to a different F\ 
width, with Dfi = 15 (top row), Dfi = 10 (middle row), and 
Dpi = 5 (bottom row). The relative exchange field orientations are 
orthogonal, with 0i = 7t/ 2, and 02 = 0. 


metal ( h/ep = 0) to a half-metallic ferromagnet ( h/ep = 1). 
It is particularly useful to consider the behavior of the sin¬ 
glet and triplet correlations over this broad range of strengths 
of h/ep. Thus, in Fig. 3, we show the same quantities as 
Fig. 2, plotted now as a function h/ep. Again, we have or¬ 
thogonal relative exchange field orientations, with 0\ = 7r/2, 
and 02 = 0. Each three-panel row corresponds to a differ¬ 
ent F\ width: Dpi = 15,10, and 5 (from top to bottom). 
The central column reveals that the averaged equal spin am¬ 
plitudes | f 1 ,avg | displays regularly occurring prominent peaks, 
the number of which varies with the length of the F\ region. 
For the exchange fields and F\ widths considered in Fig. 2, 
the triplet fi was generally weaker than either the singlet / 3 
or triplet /q. For the system parameters used in Fig. 3 how- 
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FIG. 4. (Color online) Plots of the averaged singlet and triplet com¬ 
ponents as a function of magnetic orientation /?i. Here Df i = 10, 
and results for several magnetic strengths are shown, ranging from 
weak to half-metallic. 


ever, we find that for narrow widths Dpi and sufficiently large 
exchange fields, the equal-spin triplet component /i can dom¬ 
inate the other pair correlations. In particular, for strong ferro- 
magnets with Ii/ef ~ 0.8, and thin F\ layers with Dpi = 5, 
panels (g) and (i) illustrate that the /o and fa amplitudes con¬ 
sisting of opposite spin pairs, are negligible due to the pair 
breaking effects of the strong magnet. On the other hand, the 
equal-spin pairs shown in panel (h) are seen to survive in this 
limit. This has important consequences for isolating and mea¬ 
suring this triplet component in experiments. 

Having seen how the magnitude of the exchange field h 
affects the singlet and triplet correlations, we next investi¬ 
gate the effects of changing its direction. Therefore, we ex¬ 
amine in Fig. 4, the behavior of the averaged singlet and 
triplet amplitudes when changing the magnetic orientation an¬ 
gle, /3i. We again consider a broad range of exchange field 
strengths, as shown in the legend. One of the more obvious 
features is that the maximum of |/i, a vg| typically does not 
occur for orthogonal relative exchange fields , 35 for smaller 
B\ < 90°, especially for stronger magnets. This is in agree¬ 
ment with previous 34,35,37 ’ 38 ’ 48 " 50 experimental and theoretical 
results. Due to the non-monotonicity of |/i, a vg| with h [see 
Fig. 3(e)], the h/sp = 0.35 case seen in Fig. 4(b) is larger for 
all /3i than for the weaker h/ep = 0.1 case. The singlet /3 
and triplet /o amplitudes are largest for antiparallel configu¬ 
rations (/) 1 = 180°), where the opposite exchange fields are 
effectively weakened, with reduced spin-splitting effects on 
the opposite-spin Cooper pairs. This is a well-known result. 
The results also show that the relative magnetic orientation 
angles leading to the minima of these two quantities are an¬ 
ticorrelated with the angles at which the fi correlations are 
maximal. As seen in Figs. 4(a) and (c), |/o,av g | and |/3, avg | 
decay much more abruptly as the value of h in the magnets ap¬ 
proaches the half-metallic limit: this is consistent with the dis¬ 
cussions above. Therefore, SFF structures involving strong 
ferromagnets (h ~ Ey) with f3\ at or near orthogonal ori¬ 
entations, can host larger generated triplet pair correlations 
whereby |/i| {I/0U/3I}- thus allowing for direct prob¬ 

ing of the spin triplet superconducting correlations in experi¬ 
ments. 
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FIG. 5. (Color online) Local spatial profiles of the real parts of the 
triplet components fo and /i in the F 2 region for a few different F\ 
widths, Df 1 . The exchange field in the ferromagnets corresponds to 
h/sF = 0.5, and the relative exchange field orientations are orthog¬ 
onal, with /3i = 7t/2, and f3 2 = 0. 


More detailed information regarding the triplet amplitudes, 
can be obtained from the spatial profiles of the local triplet 
correlations within the F 2 region. In Fig. 5, we present the 
real parts of the normalized fo(x) and fi(x) triplet compo¬ 
nents in terms of the dimensionless coordinate X. Results are 
plotted at four different values of Dpi as indicated in the leg¬ 
end. The exchange fields in the ferromagnets has magnitude 
corresponding to h/ep = 0.5, and the directions are mutually 
orthogonal, with j3\ = 90°, and :i 2 = 0. For the time scale 
considered here, the imaginary part of fo is typically much 
smaller than its real part. As to fi, its imaginary part is usually 
not negligible, but it exhibits trends that are similar to those 
for the real part. Examining the top panel, it is evident that 
fo exhibits the trademark damped oscillatory spatial depen¬ 
dence arising from the difference in the spin-up and spin-down 
wavevectors of the Cooper pairs. The oscillatory wavelength 
is thus governed by the quantity 27= 27r {h/ep) -1 , 
which for our parameters corresponds to 47 t. The modulat¬ 
ing fo has the same wavelength for each Dpi, although each 
curve can differ in phase. The averaged fo amplitudes are con- 
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sistent with this local behavior: Fig. 2(a) demonstrated that 
when Dpi ~ 5 and Dpi ~ 10, there is an enhancement 
of the /o component, while for Dpi « 3, it is substantially 
reduced. The equal-spin fi amplitudes, are shown in the bot¬ 
tom panel of Fig. 5. Near the interface at X = 0, the fi 
correlations are created, and then they subsequently increase 
in magnitude until deeper within the ferromagnet, where they 
clearly exhibit a gradual long-ranged decay. The trends ob¬ 
served here are opposite to those in the top panel, where for 
instance the Dpi = 3 case leads to maximal fi triplet gener¬ 
ation, in agreement with Fig. 2(b). 


Local density of states 

After the discussion of the salient features of the singlet 
and triplet pair correlations in the outer F layer, we now turn 
to the main topic of the paper: the local density of states 
measured in F 2 . This is the experimentally relevant quan¬ 
tity that can reveal the signatures of these correlations. The 
damped oscillatory behavior of the pair correlations can lead 
to spectroscopic signatures in the form of DOS inversions, 51 
and multiple oscillations. 24 In the quasiclassical approxima¬ 
tion, l4 ' 5 ' 48 - 49 a ZEP can emerge from the long-range triplet 
correlations ’ 2 -^ 2 in SFF systems. However, this approxima¬ 
tion is not appropriate for experimental conditions involving 
strong magnets and clean interfaces. It would be beneficial 
experimentally to characterize the ZEP relation to the singlet 
and triplet correlations and see how the ZEP may be a use¬ 
ful fingerprint in identifying the existence of the long-range 
triplet component. To properly do this over the broad range of 
parameters found in experimental conditions, a microscopic 
self-consistent theory that can accommodate the wide rang¬ 
ing length and energy scales is needed. In this subsection, we 
therefore present an extensive microscopic study of the ZEP 
as a function of parameters such as F layer thicknesses, ex¬ 
change energy, or interface transparency. These results are 
then correlated with the self-consistent singlet and triplet pair 
correlations in the previous subsection. In what follows, the 
DOS is normalized by the DOS at the Fermi level J\fp, and 
plotted vs the normalized energy e/Ao, where Ao is the bulk 
value of the pure S material gap at zero temperature. Our em¬ 
phasis will be on energies within the subgap region e < Ao, 
where the ZEP phenomenon arises. Since the DOS is a local 
quantity that depends on position [see Eq. (9)], in our calcula¬ 
tions we assume the location to be near the edge of the sample 
just below the STM tip as shown in Fig. 1. 

To correlate the triplet correlations in Fig. 3 with the ZEP, 
we begin by studying in Fig. 6 the sensitivity of the DOS to 
a broad range of exchange field strengths h. Each panel cor¬ 
responds to a different F\ width. Dpi. The range of h con¬ 
sidered in each panel varies since the largest ZEP depends on 
the relative values of h and Dpi. The top panel {Dpi = 5) 
clearly shows the progression of the ZEP with h: Beginning 
with the smallest exchange field, h/eF = 0.1, a moderate 
peak is observed that increases to its maximum height and nar¬ 
rower width when h/sp = 0.2. Further increases in h contin¬ 
uously diminish the ZEP, broadening its width, until eventu- 



FIG. 6. (Color online) Normalized (see text) local DOS. The multi¬ 
ple curves in each panel are for different values of H/ef- Each panel 
corresponds to a different value of Dfi (as labeled). The ferromag- 
nets have exchange fields with orthogonal relative directions. 


ally it is effectively washed away. This non-monotonic behav¬ 
ior is consistent with the ZEP being related to the presence of 
the fi triplet amplitude near the edge of the ferromagnet. This 
can be seen by reexamining the triplet amplitudes in Fig. 3(h), 
where the exchange field leading to the highest ZEP occurs 
when |/i,avg| is largest, at H/ef ss 0.2. The same consistency 
is found between Figs. 2(e) and (b) and the middle and lower 
panels of Fig. 3, respectively. For both the Dpi = 10,15 
cases, the average value of \fi\ is largest near h/ep = 0.1. 
However, the secondary peak structure in Fig. 2 is not clearly 
refected in the DOS. 

Next we study the DOS counterpart to Fig. 2. The normal¬ 
ized DOS, and the corresponding ZEP, are shown in Fig. 7 for 
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FIG. 7. (Color online) Normalized local DOS as a function of the 
normalized energy. The curves in each panel are for different val¬ 
ues of the width Dpi- Each panel corresponds to a different Ii/ef'- 
h/sF = 0.05,0.1,0.5, and we consider orthogonal relative ex¬ 
change fields. 


a broad range of widths Dpi- The parameter values here are 
similar to those used in Fig. 2, where each panel corresponds 
to a different exchange field. In panel (a) with H/ef = 0.05, 
the most prominent ZEP occurs for Dpi = 25, coinciding 
with the Fi width that yields a local maximum for the m = 0 
triplet amplitude /o (see Fig. 2(g)). By comparison, the fi 
component observed in Fig. 2(h), is smaller and lacks the mul¬ 
tiple peak structure found for /o, at this weaker exchange field. 
Therefore, the largest ZEP in the case of weak exchange fields, 
does not necessarily occur when the triplet fi is maximal; as 
Fig. 2(h) demonstrated, |/i, a vg| peaks at Dpi = 10 before 
rapidly declining. For these weaker fields, it follows from 
Fig. 2 that the magnitude of / 0 exceeds that of f 1 . It would 
appear then that it is the larger triplet component which de¬ 



FIG. 8. (Color online) Variation of the normalized local DOS with 
the in-plane exchange field angle /3i. The exchange field is fixed 
along r in F?. Also, Dfi = 10, and h/sF = 0.1. 


termines the ZEP structure. This is consistent with the known 
result 38 that the total value of the triplet component is corre¬ 
lated with T c . The next case in panel (b) corresponds also to 
a moderately weak magnet with Ii/ef = 0.1, or double the 
exchange field considered in panel (a). Since the frequency 
of the oscillations involving the opposite-spin f 0 amplitudes 
[see Fig. 2(d)] also doubles, the maximum ZEP at Dpi = 12, 
occurs at about half the h\ width found for the maximum 
ZEP in (a). The equal-spin triplet correlations fi were seen 
in Fig. 2(e) to exhibit a single peak structure, but their magni¬ 
tude is larger than at weaker fields. This is because typically 
stronger magnets in this situation lead to an enhancement of 
the fi amplitudes. Lastly, we consider (bottom panel) a rel¬ 
atively strong ferromagnet with h/ep = 0.5. For this case, 
there are additional subgap peaks flanking the main ZEP. The 
larger ZEP arises at smaller widths (Dpi = 7 . 5 , 8.5) than for 
weaker exchange fields, due to an increase in the frequency of 
the oscillations as a function of Dpi for the /o and fi com¬ 
ponents as seen in Fig. 2(a) and (b). Thus, the ZEP tends to 
exhibit a structure that dampens and widens for strong mag¬ 
nets, while the opposite is true for weaker ones and is cor¬ 
related with the stronger of the m = 0 and m = ± 1 triplet 
components present. 

Having established the behavior of the ZEP for differing 
h/ep, we now fix the magnitude of the exchange fields in 
each magnet and investigate the effects of varying their rela¬ 
tive orientation. Figure 8 illustrates the normalized DOS for 
the specific case Dpi = 10, and Ii/ef = 0.1. According 
to Fig. 4(b), the equal-spin triplet component fi is greatest 
when f3i ~ 90°. Thus we would expect the ZEP to also 
be maximal at this angle. Figure 8 shows that this is indeed 
the case. There the normalized DOS is shown for a range of 
0° < ft < 180° in increments of 30°. Clearly the orthog¬ 
onal relative exchange field (Pi = 90°) configuration results 
in the most prominent ZEP. When pi deviates from this angle 




























9 



s/Aq 


FIG. 9. (Color online) Evolution of the ZEP with scattering strength 
Hb'- The normalized local DOS is shown as a function of the di¬ 
mensionless energy. Each curve depicts results for a different scat¬ 
tering strength Hb (see text). The system parameters correspond to 
Df i = 10, and H/ef = 0.1. The exchange fields in the ferromag- 
nets are mutually orthogonal with /3i = 90°, and 02 = 0°. 


towards the P (0i = 0°) or AP (0i = 180°) alignments, both 
the triplet amplitude fi and the ZEP decline until 0i = 0° or 
180°, whereby fi = 0, and the ZEP has vanished. 

Finally, in Fig. 9 we examine the effects of interface scat¬ 
tering on the self-consistent energy spectra. We assume that 
each interface has the same delta function potential barrier 
with dimensionless scattering strength Hb- We consider a 
broad range of Hb , from transparent interfaces with Hb = 0, 
to very high interfacial scattering, with Hb = 1-6. By allow¬ 
ing Hb to vary, we effectively control the proximity effects: 
a small Hb results in stronger proximity coupling between 
the F and S regions, while a large Hb results in isolation of 
each segment, and weak proximity effects. This is evident in 
the DOS, as seen in Fig. 9, which has its largest ZEP when 
Hb = 0. The width and height of the ZEP is strongly in¬ 
fluenced by the presence of interface scattering. Increasing 
Hb results in the ZEP widening while gradually diminishing 
in height. Eventually, when the scattering strength reaches 
Hb ~ 0.7, the peak begins to split. Further increments in 
Hb causes the peaks to separate and eventually proximity ef¬ 
fects are so weakened that the DOS becomes that of an iso¬ 
lated bulk ferromagnet. The two secondary subgap peaks that 
lie symmetrically about the ZEP are seen to also decline in a 
monotonic fashion as Hb becomes larger. 


B. Diffusive Regime 

In this section, we consider a diffusive SFF junction in 
the full proximity limit. We employ the Usadel approach de¬ 
scribed in Sec. II to investigate the local DOS. As remarked 


earlier, the quasiclassical method is limited to energies close 
to the Fermi level. Hence, our discussion here will be limited 
to relatively weak ferromagnets. As in the ballistic regime, 
we consider heterostructures where the magnetic layers are 
made of identical materials so that the ferromagnetic coher¬ 
ence lengths are the same, = £p 2 = t;p, and we consider 
the low temperature regime where T = 0.05T C . Prior to cal¬ 
culating the DOS, we normalize the Usadel equation by 
which in the diffusive regime is written, = yjD/h. Us¬ 
ing this normalization scheme, the explicit dependency on the 
exchange field is removed and the Usadel equation now in¬ 
volves terms containing the ratio dp/^p. This approach can 
lead to easier pinpointing of regions in parameter space where 
the ZEP is most prominent, and it also permits a broad range 
of this ratio to be studied. We assume that the magnetic orien¬ 
tation angle is fixed at /3 2 = 0, or equivalently h = (0,0 ,h z ). 

We numerically solve the Usadel equation, Eq. (10), to¬ 
gether with the mentioned boundary conditions. To find the 
total Green’s function, we substitute the solution into Eq. (10) 
and obtain the DOS. To determine the optimal geometry in 
which the ZEPs are most pronounced, we present in Fig. 10 
the ZEP at the topmost edge of the SFF structure, corre¬ 
sponding to the location x = dpi -\-dp 2 - The two-dimensional 
color mapping depicts the strength of the DOS at zero en¬ 
ergy (the ZEP) as a function of the normalized F thicknesses, 
dp \ //f and dP 2 /0F- In the top row panels, the internal field 
of the F layers have a misalignment angle of B-\ = 7r/2, 
while for the bottom row 0i = tt/ 6. The left, middle, and 
right columns are for different opacities at the SF interface: 
C, = 1.0, 2.5, and 5.0, respectively. By increasing Q the over¬ 
all strength of the proximity effects is effectively weakened: 
it is evident that transparent SF contacts yield stronger ZEPs, 
that persist in thicker F layers. It is also apparent that the 
orthogonal case 0i = 7 t /2 has more extensive regions in the 
parameter space spanned by the F thicknesses with enhanced 
ZEPs, as compared to the fii = n/6 case. An important as¬ 
pect of the ZEP that all cases investigated in Fig. 10 share, is 
that it is strongest when dpi <C dp 2 - This finding is fully con¬ 
sistent with low proximity bilayer SFF hybrids. 48 Therefore, 
for the parameters considered here, the ZEPs are strongest for 
C, = 1, 0.5£f S dpi < £p , and 1.5 £p 0/ dp 2 3.5^. 
The ratio of the F thickness to the length scale is an im¬ 
portant dimensionless quantity that appears in the normalized 
Usadel equations, and consequently thinner dpi and dp 2 , al¬ 
low for stronger ferromagnets when studying the DOS in the 
diffusive limit. 

Finally, we study the sensitivity of the ZEPs to both the ori¬ 
entation angle 0i and interface transparency parameter (/. We 
thus show in Fig. 11 the ZEP as a function of 0i over a wide 
range of /, as shown in the legend. The geometric parameters 
correspond to dpi = 0.85£f and dp 2 = 3.5£f, which resides 
within the range of system widths studied in Fig. 10 resulting 
in the largest ZEPs. In calculating the ZEP, we again consider 
the DOS at the edge of the sample (see also Fig. 1). It is seen 
that the interface transparency can significantly alter the be¬ 
havior of the ZEP as the relative exchange field angle sweeps 
from the P (/3 2 = 0°) to AP (/? 2 = 180°) orientations. For 
example, when (' = 1.0, the maximal ZEP is offset from the 
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FIG. 10. (Color online) Zero energy peak 
in the DOS spectrum of diffusive SFF spin 
valves as a function of the normalized F layer 
thicknesses d F i/£ F and d F 2 /(,!•'■ Each col¬ 
umn corresponds to a different SF trans¬ 
parency, with £ = 1,2.5,5 (from left to 
right). The top row of panels shows the evo¬ 
lution of the ZEP for the misalignment angle 
Pi = 7t/2, while the bottom row of panels 
are for Pi = n/Q. For both cases, the in¬ 
ternal field of the F2 layer is along the z di¬ 
rection, P 2 = 0. The ZEPs are computed at 
x = d F i + d F 2 (at the top most F /vacuum 
interface). 



FIG. 11. (Color online) Zero energy peak of the DOS spectrum in 
diffusive SFF spin valves as a function of exchange field orientation 
pi for several values of £, which controls the opacity of the SF 
interface. We set P 2 = 0 and rotate the exchange field direction of 
F\ from the parallel (Pi = 0) to antiparallel (Pi = ir) orientations. 
We have chosen representative values of d F 1 = 0.8£f and d F 1 = 
3.5£f, in accordance with the system parameters used in Fig. 10. 


orthogonal configuration, occurring at P 2 ~ 0.67T. By increas¬ 
ing the barrier strength, this peak shifts towards larger P 2 , un¬ 
til the relative exchange fields are nearly antiparallel. There 
is also a simultaneous reduction in amplitude, due to the F 
and S regions becoming decoupled as the proximity effects 
diminish. Interestingly, as ( increases, there is a splitting of 
the main peak: weaker secondary peaks emerge. Eventually 
however, for sufficiently large £, the opacity of the interface 
causes the low energy DOS to be insensitive to p 1 , and the 
ZEP flattens out. The ZEPs are also observed to disappear 
when the relative exchange fields are collinear, corresponding 
to the situation when the triplet amplitudes vanish in both the 
diffusive and ballistic regimes (see also Fig. 4). 


IV. SUMMARY AND CONCLUSIONS 

In summary, we have employed a microscopic self- 
consistent wavefunction approach to study the low energy 
proximity induced local DOS in SFF spin valves with non- 
collinear exchange fields in the clean limit. Our emphasis has 
been on the results of STM methods that probe the outer F 
layer. To identify the physical source of the corresponding 
ZEPs that occurs for such data in these systems, we also cal¬ 
culated the absolute value of the triplet pair correlations, av¬ 
eraged over the outer F layer. We have done so for a broad 
range of experimentally relevant parameters, including the ex¬ 
change field strength and orientation, as well as thicknesses 
of the ferromagnets. Our results demonstrate a direct link 
between the spin-1 triplet correlations and the appearance of 
ZEPs in the local DOS spectra, and point to system parame¬ 
ters and configurations which would support larger equal-spin 
triplet superconducting correlations. These correlations could 
then be probed indirectly via single-particle signatures that are 
measurable using local spectroscopy techniques. Our results 
are consistent with’ 8 findings relating the average strength of 
triplet correlations to the angular dependence of the transi¬ 
tion temperature. Our findings suggest that the ZEPs arising 
from the spin-1 triplet amplitudes can be effectively isolated 
in SFF systems with strong ferromagntets, with the outer one 
being very thin. This asymmetric geometry not only produces 
greater equal-spin triplet generation, but it can also filter out 
the rapidly decaying opposite-spin pairs deep within the sam¬ 
ple. We also considered the same valve structure in the diffu¬ 
sive regime utilizing a Green function method within the full 
proximity limit. Our investigations yielded a broad range of 
F layer thicknesses and relative exchange fields orientations 
that lead to observable signatures in the low energy DOS, thus 
also giving useful guidelines for future experiments. 
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